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Abstract 

Macroscopic traffic models liave recently been severely criticized to base on 
lax analogies only and to have a number of deficiencies. Therefore, this paper 
shows how to construct a logically consistent fluid-dynamic traffic model from 
basic laws for the acceleration and interaction of vehicles. These considera- 
tions lead to the gas-kinetic traffic equation of Paveri-Fontana. Its stationary 
and spatially homogeneous solution implies equilibrium relations for the 'fun- 
damental diagram', the variance-density relation, and other quantities which 
are partly difficult to determine empirically. 

Paveri-Fontana's traffic equation allows the derivation of macroscopic mo- 
ment equations which build a system of non-closed equations. This system 
can be closed by the well proved method of Chapman and Enskog which leads 
to Euler-like traffic equations in zeroth-order approximation and to Navier- 
Stokes-like traffic equations in first-order approximation. The latter are finally 
corrected for the finite space requirements of vehicles. It is shown that the 
resulting model is able to withstand the above mentioned criticism. 
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I. INTRODUCTION 



Because of analogies with gas theory |]T]-0 and fluid dynamics P-p|,|3|,p!0| modeling and 
simulating traffic flow increasingly attracts the attention of physicists , 11 - 14 ] . How- 

ever, due to the great importance of efficient traffic for modern industrialized countries, 
the investigation of traffic flow has already a long tradition. In the 1950s Lighthill and 
Whitham |TD| as well as Richards [|I^] proposed a first fluid- dynamic (macroscopic) traffic 
model. During the 1960s traffic research focused on microscopic follow-the-leader models 
||TB|^H[. Mesoscopic models of a gas-kinetic (Boltzmann-like) type came up in the 1970s 
2^ , |25| flj3|J^,P^ . Since the 1980s simulation models [^,^] play the most important role due 



to the availability of cheap, fast, and powerful computers. We can distinguish macroscopic 



traffic simulation models P9H3^, microscopic simulation models ||33H36| which include cel- 



lular automaton models p7 H39|JTT| -|l4| , and mixtures of both 

In high-fidelity microscopic traffic models each car is described by its own equation(s) 
of motion. Consequently, computer time and memory requirements of corresponding traffic 
simulations grow proportional to the number of simulated cars. Therefore, this kind 
of models is mainly suitable for off-line traffic simulations, detail studies (for example of 
on-ramps or lane mergings), or the numerical evaluation of collective quantities like 
the density-dependent velocity distribution, the distribution of headway distances etc., and 
other quantities that are difficult to determine empirically. 

For this reason, fast low-fidelity microsimulation models that allow bit-handling have 
been developed for the simulation of large freeways or freeway networks []37| , ^ . However, 
although they reproduce the main effects of traffic flow, they are not very suitable for detailed 
predictions because of their coarse-grained description. 

Therefore, some authors prefer macroscopic traffic models |ll0tp^ - |43|j30| , |4^ j5|-|9| . These 
base on equations for collective quantities like the average spatial density p{r, t) per lane (at 
place r and time t), the average velocity K(r, t), and maybe also the velocity variance G(r, t). 
Here, simulation time and memory requirements mainly depend on the discretization Ar 
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and At of space r and time t, but not on the number of cars. Therefore, macroscopic 
traffic models are suitable for real-time traffic simulations. The quality and reliability of the 
simulation results mainly depend on the correctness of the applied macroscopic equations 
and the choice of a suitable numerical integration method. The rather old and still continuing 
controversy on these problems pT] , ^^^^^pD| , P| , ^5| -pIl|,|8|,|9|,|5|JI|J^ shows that they are 



not at all trivial. 

Some of the most important points of this controversy will be outlined in Section II. It 
will be shown that even the most advanced models still have some serious shortcomings. 
The main reason for this is that the proposed macroscopic traffic equations were founded on 
heuristic arguments or based on analogies with the equations for ordinary fluids. In contrast 
to these approaches, this paper will present a mathematical derivation of macroscopic traffic 
equations starting from the gas-kinetic traffic equation of Paveri-Fontana which is very 



reasonable and seems to be superior to the one of Prigogine and co-workers . The 

applied method is analogous to the derivation of the Navier-Stokes equations for ordinary 
fluids from the Boltzmann equation It bases on a Chapman-Enskog expansion | pB|JF7| 

which is known from kinetic gas theory and leads to idealized, Euler-like equations in zeroth- 
order approximation and to Navier-Stokes-like equations in flrst-order approximation 



In this respect, the paper puts into effect the method suggested by Nelson A similar 
method was already applied to the derivation of fluid-dynamic equations for the motion of 
pedestrian crowds but it assumed some dissatisfactory approximations. 



The further procedure of this paper is as follows: Section II presents a short history of 
macroscopic traffic models and discusses the abilities and weaknesses of the different ap- 
proaches. Section III introduces the Boltzmann-like model of Prigogine |^ and compares 
it with the one of Paveri-Fontana 0. From their gas-kinetic equations macroscopic ('fluid- 
dynamic') traffic equations will be derived in Section IV. Unfortunately, they turn out to 
build a hierarchy of non- closed equations, i.e. the density equation depends on average veloc- 
ity V, the velocity equation on velocity variance G, etc. Therefore, a suitable approximation 
must be found to obtain a set of closed equations. It will be shown that some of the traffic 



models introduced in section II correspond to zeroth-order approximations of different kinds. 
These, however, are not very well justified. A similar thing holds for the Euler-like traffic 
equations which, apart from a complementary covariance equation, contain additional terms 



compared with the Euler equations of ordinary fluids ||58|. These are, on the one hand, due 
to a relaxation term which describes the drivers' acceleration towards their desired velocities. 
On the other hand, they are due to interactions which are connected with deceleration pro- 
cesses since these do not satisfy momentum and energy conservation in contrast to atomic 
collisions. 

A very realistic, first-order approximation which is, in a certain sense, self-consistent can 
be found by solving the reduced Paveri-Fontana equation which is obtained from the original 
one by integration with respect to desired velocity. We will utilize the fact that, according 
to empirical traffic data |]60| , |6l| ,|3|, |62| , |33| , the equilibrium velocity distribution has a Gaussian 



form. This allows the derivation of mathematical expressions for the equilibrium velocity- 
density relation, the 'fundamental diagram' of traffic flow, and the equilibrium variance- 
density relation (cf. Sec. IV C). Afterwards an approximate time-dependent solution of 
Paveri-Fontana's equation will be calculated by use of the Euler-like equations. Due to 
the additional terms in Paveri-Fontana's equation compared with the Boltzmann equation 
the corresponding mathematical procedure is more complicated than the Chapman- Enskog 
expansion for ordinary gases (cf. Sec. V). 

Nevertheless, it is still possible to derive correction terms of the Euler-like macroscopic 
traffic equations (cf. Sec. VI). These have the meaning of transport terms (like e.g. the 
flux density of velocity variance) and are related with the finite skewness 7 of the velocity 
distribution in non-equilibrium situations. The resulting equations are Navier-Stokes-like 
traffic equations which, in comparison with the ordinary Navier-Stokes equations [0], con- 
tain additional terms arising from the acceleration and interaction of vehicles. Additionally, 
they are complemented by a covariance equation which takes into account the tendency of 
drivers to adapt to their desired velocities. 

Because of the one-dimensionality of the Navier-Stokes-like traffic equations no shear 



viscosity term occurs. However, in Section VII it is indicated how transitions between 
different driving modes can cause a bulk viscosity term. Furthermore, corrections due to 
finite space requirements of each vehicle (vehicle length plus safe distance) are introduced. 

The resulting model overcomes the shortcomings of the former macroscopic traffic models 
(that are mentioned in Sec. II). Section VIII summarizes the results of the paper and gives 
a short outlook. 



II. SHORT HISTORY OF MACROSCOPIC TRAFFIC MODELS 



In 1955 Lighthill and Whitham proposed the first macroscopic (fluid-dynamic) traffic 
model. This bases on the continuity equation 

dp d{pV) 



dr 







(1) 



which reflects a conservation of the number of vehicles. For the average velocity V , Lighthill 
and Whitham assumed a static velocity-density relation: 



V{T,t) := V;[p(r,t)]. 



(2) 



Inserting (^ into ([T|) we obtain 



dp 



Ve + P 



dp 



or 



(3) 



Equation (^ describes the propagation of non-linear 'kinematic waves ' with velocity c(p) = 
Ve{p) + pdVe/dp p]| , |55| . In the course of time the waves develop a shock structure, i.e. their 
back becomes steeper and steeper until it becomes perpendicular, leading to discontinuous 
wave profiles ||IO|Jl5| , |63[ . 



In reality, density changes are not so extreme. Therefore, it was suggested to add a 
diffusion term Dd'^p/dr'^ which smoothes out the shock structures somewhat |]63| , |6^ . The 
resulting equation reads 

dp 



dp 
dt 



^ dr 



dVe dp ^d'^p 



dp dr 
5 



Q^2 



(4) 



For the case of a linear velocity-density relation |65 



K(P) :=Knaxfl-^) (5) 

\ Pmax / 



it can be transformed into the Burgers equation 

dg dg £)^^9 (q\ 
dt ^ dr dr"^ 

which is analytically solvable [^. Here, we have introduced the function 



^7[p(r,t)]:=Knaxfl-^^) . (7) 

\ Pmax / 

The most important restriction of models (|l]), @) and @, (H) is relation which 
assumes that average speed V{r,t) is always in equilibrium with density p{r,t). Therefore, 
these models are not suitable for the description of non-equilibrium situations occuring at 
on-ramps, changes of the number of lanes, or stop-and-go traffic. 

Consequently, it was suggested to replace relation (0) by a dynamic equation for the 
average velocity V. In 1971, Payne ||5l[] introduced the velocity equation 

dV ^BV C(p)dp K,,, , 

with 



C(p). ''''' ' 



dp 



(8b) 



2r dp 2t 

which he motivated by a heuristic derivation from a microscopic follow-the-leader model 
|67| . Here, VdV/dr is called the 'convection term' and describes velocity changes at place 
r that are caused by average vehicle motion. The 'anticipation term' —{C/p)dp/dr was 
intended to account for the drivers' awareness of the traffic conditions ahead. Finally, the 
'relaxation term' [Ve(p) — V]/t delineates an (exponential) adaptation of average velocity V 
to the equilibrium velocity Ve{p) with a relaxation time t. 

Unfortunately, for bottlenecks the corresponding computer simulation program 'FRE- 
FLO' suggested by Payne produces output that "does not seem to reflect what really 
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happens even in a qualitative manner" As a consequence, several authors have sug- 

gested a considerable number of modifications of Payne's numerical integration method or 
of his equations p8| , ^ , ^ , ^ , |30| , ^ , ^ , ^ , |69| . A more principal weakness of Payne's equa- 
tions is that their stationary and homogeneous solution is stable with respect to fluctuations 
over the whole density range which can be shown by a linear stability analysis p5| , ^ , P 



Therefore, Payne's model (|T|), (H) does not describe the well-known self-organization of stop- 
and-go waves above a critical density |^J7^. This problem is removed by substituting 
relation (|bD by 



C{p) 



dp 



(9) 



with the equilibrium 'traffic pressure ' 



Ve{p) := pQe{p) . 



(10) 



The modifled velocity equation reads 



dt Or 



:ila) 



p dr T 

and can be derived from the gas-kinetic (Boltzmann-like) traffic models 0,0,0 (cf. Section 
IV). For 6e(p), Phillips |3|,[7T1| suggested a relation of the form 



<de{p) := Oo 1 



P 



(lib) 



In contrast, Kiihne as well as Kerner and Konhauser p],^ assumed, as a first approach, 
6e to be a positive constant: 



Qeip) ■■= eo. 



(12) 



Unfortunately, equations (Q), ( |lla] ) predict the formation of shock waves like Lighthill and 
Whitham's equation does P3|,p|. For this reason, Kiihne ||i3| , |70[] suggested to add a small 
viscosity term ud'^V/dr'^ which smoothes out sudden density and velocity changes somewhat. 
Then, the velocity equation 



^^V—-- — — + u— ' ^ 
dt dr p dr dr 



+ f^ = --^ + ''^ + ;IW-fl (13) 



results. A linear stability analysis of Kiihne's equations (|I]), (|l^) shows that these predict 
the self-organization of stop-and-go waves or of so-called 'phantom traffic jams' (i.e. unstable 
traffic) on the condition 



Pe 



> a/0o(1 + rz/A;2) (14) 



dp 

where k denotes the wave number of the perturbation |]73|,|5| . This condition is fulfilled if the 



equilibrium density pe corresponding to the stationary and spatially homogeneous solution 
exceeds a critical density per that depends on the concrete form of Ve{p). 

For reasons of compatibility with the Navier-Stokes equations for ordinary fiuids Kerner 
and Konhauser replaced Kiihne's constant u by the density- dependent relation 

v{p) = J (15) 

with the constant viscosity coefficient uq. Computer simulations of their equations (|l]) and 
(p!3|), (|15D show the development of density clusters if the critical density pcr given by 



(|l^) and ([TsD is exceeded. On the basis of a very comprehensive study of cluster-formation 
phenomena, Kerner and Konhauser presented a detailed interpretation of stop-and-go 
traffic. 

Despite the considerable variety of proposed macroscopic traffic models, even the most 
advanced of them have still some shortcomings. For example, for a certain set of parameters 
the mentioned models predict traffic densities that exceed the maximum admissible density 
Pbb = 1/^0 which is the bumper-to- bumper density {Iq = average vehicle length) ||5|. Fur- 
thermore, in certain situations even negative velocities may occur To illustrate this, 



imagine a queue of vehicles of constant density po- Assume that, e.g. due to an accident 
that blocks the road, this queue has come to rest (i.e. ^ = 0) and that it ends at r = tq 
which shall imply p(r, t) = for r < tq. Then, dp/dr diverges at place tq (or is at least very 
large) and equations (D, (lUD, (|13D all predict dV{ro,t)/dt < if 6 ^ 0. 
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Of course, we wish to have a model that is not only valid in standard situations, but also 
in extreme ones. Moreover, the model should provide reasonable results not only for certain 
parameter values. This is particularly important for the reason that technical measures like 
automatic distance control may change some parameter values considerably. Nobody knows 
if the existing phenomenological models are still applicable, then. Therefore, we will derive 
the specific structure of the traffic model from basic principles regarding the behavior of the 
single driver-vehicle units and their interactions. 

III. GAS-KINETIC (BOLTZMANN-LIKE) TRAFFIC MODELS 

Let us assume that the motion of an individual vehicle a can be described by several 
variables like its place rQ,(t), its velocity fQ,(t), and maybe other quantities which characterize 
the vehicle type or driving style (the driver's personality). We can combine these quantities 
in a vector 



that denotes the state of vehicle a at a given time t. The time- dependent phase-space density 



is then determined by the mean number An{r,v, . . . ,t') of vehicles that are at a place 
between r — Ar/2 and r + Ar/2, driving with a velocity between v — Av/2 and v + Av/2, 
. . . at a time t' E [t - At/2, t + At/2]: 

t+At/2 

1 r 

p{r,v,...,t)ArAv...:=— J dt' An{r,v, . . . ,t') . (18) 

t-At/2 

For vehicles, the phase-space densitiy /) is a very small quantity. Therefore, in the limit 
Ar 0, Av — > 0, . . ., At ^ it is only meaningful in the sense of the expected value of 
an ensemble of macroscopically identical systems The interpretation of p as a quantity 
which can describe single traffic situations is only possible for "coarse-grained averaging" 




(16) 



p(f,t) = p{r,v, ...,t) 



(17) 
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where Ar, Av, . . ., and At must be chosen "microscopically large but macroscopically small" 
p|.|59[] or, more exactly, 



1. smaller than the scale on which variations of the corresponding macroscopic quantities 
occur, 

2. so large that An{r, v, . . . ,t) ^ 1 which is not always compatible with the first condi- 
tion. 

However, in any case a suitable gas-kinetic equation for the phase-space density p allows the 
derivation of meaningful equations for collective ('macroscopic') quantities like the spatial 
density p{r,t) per lane, the average velocity V{r,t), and the velocity variance 6(r, t). To 
obtain an equation of this kind, we will bring in the well-known fact that the temporal 



evolution of phase-space density p is given by the continuity equation 

which again describes a conservation of the number of vehicles, but this time in phase-space 
Q = {all admissible states x}. Whereas (pdx/dt) reflects changes of phase-space density 
p due to a motion in phase space Q with velocity dx/dt, the term (9p/9t)tr delineates changes 
of p due to discontinuous transitions between states. 

A. Prigogine's model 

In Prigogine's model the state x is given by the place r and velocity v = dr/dt of a. vehicle. 
The transition term {dp/dt)tT consists of a relaxation term {dp/dt)j.ei and an interaction term 



{dp/dt)int [H^il^il- Therefore, equation ( pj] ) assumes the explicit form 



1-^^ I; (^a- (!),/©,„.■ (-) 

The interaction term {dp/dt)\nt is intended to describe the deceleration of vehicles to the 
velocity of the next car ahead in situations when this moves slower and cannot be overtaken. 
Prigogine l2^,§] suggests to describe processes of this kind by the Boltzmann equation 



10 



■.= Jdw{l-p)\v-w\p{r,v,t)p{r,w,t) (21a) 

' int y 

V 

dw (1 — p)\w — v\p(r, w, t)p(r, v, t) (21b) 

oo 

= (1 — J9)p(r, f , t) j dw {w — v)p{r,w,t) . 



where p denotes the probabihty that a slower car can be overtaken. Functional relations for 

p^pip,V,e) (22) 
are proposed in Refs. Pfl|75| . The term ( pia| ) corresponds to situations where a vehicle 



with speed w > v must decelerate to speed v, causing an increase of phase-space density 
p(r, V, t). The rate of these situations is proportional 

1. to the probability {1 — p) that passing is not possible (which corresponds to the 'scat- 
tering cross section' in kinetic gas theory), 

2. to the relative velocity |f — w| of the interacting vehicles, 

3. to the phase-space density p(r, v, t) of vehicles which may hinder a vehicle with velocity 

w > V, and 

4. to the phase-space density p{r,w,t) of vehicles with velocity w > v that may be 
affected by an interaction. 



Term ( pib|) describes a decrease of phase-space density p(r, t>, t) due to situations in which 



vehicles with velocity v must decelerate to a velocity w < v. A more detailed discussion of 
interaction term (|21| ) can be found in Refs. 

Note that approach ( |2lD assumes an instantaneous adaptation of velocity which does not 
take any braking time. Moreover, the deceleration process of the faster vehicle is assumed 
to happen at the location r of the slower vehicle, i.e. vehicles are implicitly modelled as 
point-like objects without any space requirements. The first assumption is only justified 
for braking times that are short compared to temporal changes of phase-space density p, 
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but modifications for finite braking times are possible Tlie second assumption is only 



acceptable for very small densities at which the average headway distance is much larger 
than average vehicle length plus safe distance. It will, therefore, be corrected in Section 
VII. The corresponding modifications also implicitly take into account the pair correlations 
of succeeding vehicles |]7B|. These are neglected by approach ( PT| ) due to its assumption of 



'vehicular chaos according to which the velocities of vehicles are not correlated until they 
interact with each other 

Now, we come to the description of acceleration processes by vehicles that do not move 
with their desired speeds. In this connection, Prigogine suggests a collective relaxation of 
the actual velocity distribution 

n.r.y.-^ (.3) 

towards an equilibrium velocity distribution Po{v) instead of an individual speed adjustment 
so that 

dv , 

In detail, Prigogine starts from the observation that free traffic is characterized by a certain 
velocity distribution Po{v) which corresponds to the distribution Po{vo) of desired velocities 
Vq. Moreover, he assumes that the drivers' intention to get ahead with their desired speeds 
causes the phase-space density p(r, v, t) to approach the equilibrium phase-space density 

Mr,v,t) := pir,t)Po{v) (25) 

(exponentially) with a certain relaxation time r which is given by the average duration of 
acceleration processes. Therefore, Prigogine's relaxation term has the form [p^,p^,^ 

{dp\ p{r,t)Po{v) - p{r,v,t) 



dt 



(26) 



rel 

Despite the merits of Prigogine's stimulating model, this approach has been severely crit- 
icized p|,|5l||. In a clear and detailed paper P| Paveri-Fontana showed that Prigogine's model 
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has a number of peculiar properties which are not compatible with empirical findings. For 
example, he demonstrates that the relaxation term (^) corresponds to discontinuous veloc- 
ity changes which take place with a certain, time-dependent rate. Furthermore, Daganzo 
criticized that, according to (p6D, "the desired speed distribution is a property of the road 
and not the drivers" |^ which was already noted by Paveri-Fontana [Q. In reality, however. 



one can distinguish different 'personalities' of drivers: 'aggressive' ones desire to drive faster, 
'timid' ones slower. Therefore, Paveri-Fontana developed an improved gas-kinetic traffic 
model which corrects the deficiencies of Prigogine's approach. 

B. Paveri-Fontana's model 

Paveri-Fontana assumes that each driver has an individual, characteristic desired velocity 
vq- Consequently, the associated states x are given by place r, velocity v, and desired 
velocity vq so that Prigogine's phase-space density p{r,v,t) is replaced by p{r,v,vo,t). The 
corresponding gas-kinetic equation (|19]) explicitly reads [ [T^ 

dp dijw) d ( Jv\ d ( Jvo\ (dp\ 

The term d{pdvQ/dt)/dvQ can be neglected since the desired velocity of each driver is nor- 
mally time-independent during a trip which implies 

In contrast to Prigogine, Paveri-Fontana describes the acceleration towards the desired ve- 
locity vo by 

-:=-(„o-.) (27c) 



which means an individual instead of a collective relaxation. Relation ( |27c| ) can be easily re- 
placed by other acceleration laws dv/dt or density-dependent driving programs as suggested 
by Alberti and Belli |2^. Alternatively, for acceleration processes an interaction approach 



can be formulated which was recently proposed by Nelson |T[]. However, the assumption 
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( P7c| ) of exponential relaxation is a relatively good approximation since drivers gradually 
reduce the acceleration as they approach their desired velocity fo- 

Paveri-Fontana needs the transition term (i9p/i9t)tr only for the description of deceleration 
processes due to vehicular interactions. For these he assumes the Boltzmann equation 

i^-p) J dw J dwo \v - w\p{r,v,Wo,t)p{r,w,Vo,t) 

V / tr y 

V 

— (1— p) / dw / dwo\w — v\p{r,w,wo,t)p{r,v,vo,t) (27d) 







which has an analogous interpretation as (^Tp. (For details cf. Ref. 0.) Note that, according 
to (|27dD , "the velocity of the slow car is unaffected by the interaction or by the fact of being 
passed" and that "no driver changes his desired speed" during interactions. Therefore, 



the interaction term (|27d|) fulfils the requirements called for by Daganzo ||51| : 



1. that "a car is an anisotropic particle that mostly responds to frontal stimuli" |EI 



and that "a slow car should be virtually unaffected by its interaction with faster cars 
passing it (or queueing behind it)" |]5T| . 



2. that "interactions do not change the 'personality' (aggressive/timid) of any car" [^. 



Finally, note that the proportion of vehicles jamming behind slower cars cannot accelerate. 
This circumstance can be taken into account by a density- and maybe velocity- or variance- 
dependence of the relaxation time B,|3|j75[|: 



r^rip,V,e). (28) 

In order to compare Paveri- Font ana's traffic equation with Prigogine's one we integrate 
equation (|27|) with respect to Vq and obtain the reduced Paveri-Fontana equation 



dp ^ d{vp) ^ d 
dt dr dv 



p{r,v,t) 



{1 — p)p{r,v,t) / dw {w — v)p{r,w,t) . (29) 







Here, we have introduced the reduced phase-space density 
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p(r, V, t) := J dvo p{r, v, Vq, t) (30) 

and the quantity 

Voiv;r,t):=jdvoVo-^^. (31) 

The only difference with respect to Prigogine's formulation (pOD to (^) is obviously the 
other relaxation term. 

IV. DERIVATION OF MACROSCOPIC TRAFFIC EQUATIONS 

Since we are mainly interested in the temporal evolution of collective ('macroscopic') 
quantities like the spatial density 



p{r,t) := J dvp{r,v,t) (32) 

per lane, the average velocity 

V(,t).W:=/d„.^. (33) 

and the velocity variance 

e(r, t) = {[v- V{r, t)]'') := I dv[v- V{r, t)] 



2\ . /^..L. .M2pir,v,t) 



p{r,t) 

= {v')~[Vir,t)r (34) 
we will now derive equations for the moments rrikfi with 

mk,i{r,t) = p{r,t){v^{vo)^) ■= dv dvov'^ivoY p{r,v,Vo,t) . (35) 



By multiplying Paveri-Fontana's equation ( pQ]) with and integrating with respect to v we 
obtain via partial integration, 



d d f , k 9 f J 

-Ql^kfi + + j dvv ^ I P- 

-Ql^kfl + - J dvkv I p 
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d d k 

w:mkfl + 7r"^fc+i,o {mk-i,! - nikfi) (36a) 

ot or T 



(1 — p) J dvp{r,v,t) J dw (wv^ — v'^^^)p{r,w,t) 

(1 - p) (mi,omfc,o - ?^fc+i,o"^o,o) • (36b) 



Applying the analogous procedure to Prigogine's model ( poD to (|26D , for the moments 



rrikflir, t) = p{r, t) {v^) := j dv v^p{r, v, t) (37) 



one can derive the equations 



d d 1. 

TTi^kfi + 7r"^fe+i,o = -{iriQ^k - rrikfi) 
ot or T 

+ (1 - p)(mi,omfc,o - "^fe+i,o"^o,o) (38) 



(cf. [H) where 



mo^k{r,t) := J dvo {vo)''po{r,Vo,t) 

= p(r,t) I dvoivo)'Poivo). (39) 

A comparison of moment equations (|36D with (|38[) shows that Prigogine's and Paveri- 
Fontana's model lead to identical equations for spatial density p(r, t) = mofi{r, t) and average 
velocity V{r,t) = mi^Q{r,t)/ p{r,t), despite the different approaches for the relaxation term. 
However, the equations for higher order moments mkfi{r,t) with k>2 differ. 

Obviously, equations (|36D as well as (|3^) represent a hierarchy of non-closed equations 
since the equation for the kth moment q depends on the {k + l)st moment mfc+i q. As a 
consequence, the density equation 

depends on average velocity V, the velocity equation 

f-f-^^-Vo-.)-a-..e 
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on variance 0, etc. Here, we have introduced the average desired velocity 

p{r,v,Vo,t) 



Vo{r,t) := dv dvoVo 



p{r,t) 



(42) 



the so-called 'traffic pressure ' P3|,|3|J7I 

1 



V{r,t) :-- 



dv {v — V)p{r, V, t) / dw {v — w)p{r, w, t) 



p{r,t) 

dv {v - Vypir, V, t) = p{r, t)Q{r, t) 



(43) 



and the equilibrium velocity 

K(p, V, 9) := Vo - T{p, V, 6) [1 - p{p, V, Q)]V 

which is related with stationary and spatially homogeneous traffic flow. 

Equations ( ^O]) and (|4ll) are easily derivable from the moment equations 
respectively by use of mo,o = P and 

dmi^o d{pV) dV dp 
dt ^ dt ^^~dt^ dt' 



(44) 



and fpSf) 



(45) 



The variance equation is obtained analogously. For the traffic equation of Paveri-Fontana it 
reads 



dQ ^dQ 
dt dr 



^^dV IdjpT) 
dr p dr 

+ -{c-e)-{i-p)pT 

T 

2V dV 1 dj 

p dr p dr 

+ -[ee(p,i^,e,c,:r)-0] 

r 



(46) 



and depends on the covariance 



C{r,t) = {{v-V){v,-V,)) 

dvo J dv{v - V){vq - Vq) 

dv {v-V)[Vo{v)-Vo 



p{r,v,Vo,t) 
p{r,t) 
p{r,v,t) 
p{r,t) 



(47) 
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r(r, t) ^ {{v - VY) := Jdv{v- v f-^^^ ■ (48) 
In addition, we have introduced the flux density of velocity variance 

J{r, t) := — [ dv{v — V)^p{r, v, t) f dw (f — w)p{r, w, t) 



as well as the third central moment 



p{r,t) 

dv {v - Vfp{r, V, t) = p{r, t)T{r, t) (49) 



(which corresponds to the 'heat flow' in conventional fluid-dynamics) and the equilhrium 
variance 

e,(p, V, e, c, j):=c- <p2Xi^ii _ ^(p, e)] j . (so) 



A. Approximate closed macroscopic traffic equations 

We will now face the problem of closing the hierarchy of moment equations by a suitable 
approximation. The simplest approximations replace a macroscopic quantity Q{r,t) (which 
would be determined by a dynamic equation) by its equilibrium value Q^. which belongs to 
the stationary and spatially homogeneous solution. Approximations of this kind are zeroth- 
order approximations. The simplest one is obtained by a substitution of V{r, t) (which 
actually obeys Eq. (^)) by the equilibrium velocity 

Ve{p) ■■= Vo - r(p)[l -p(p)]pe,(p) (51) 

(cf. (0)). Equations (^Oj), (|5TD obviously correspond to the model (|lD, (||) of Lighthill and 
Whitham. Relation ( pTD specifles the equilibrium velocity-density relation (0) in accordance 
with Paveri-Fontana's traffic equation. It could be interpreted as a theoretical result con- 
cerning the dependence of Ve{p) on the microscopic processes of traffic flow: According to 
(pl|), the equilibrium velocity K(p) is given by the average desired velocity Vq diminished 
by a term arising from necessary deceleration maneuvers due to interactions of vehicles. 
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However, according to equation ( pT]) , the approximation V{r,t) ^ V"e[p(r, t)] is only 
justified for r ^ wliicli is not compatible witli empirical data. Consequently, tlie latter 
does not adequately describe non-equilibrium situations like on-ramp traffic or stop-and-go 
traffic where the velocity is not uniquely given by the spatial density p{r,t). 

Another zeroth-order approximation is found by leaving Eq. (^TJ) unchanged but replac- 
ing the dynamic variance 0(r, t) by the equilibrium variance 

Qeip, V) := C^ip, V) - ^^^[1 - p{p, V)]pT,{p, V) (52) 

(cf. ([50|)). (Here, the subscript e shall again indicate the equilibrium- value or -relation of a 
function.) The resulting model (POD, (|^, (p^) obviously corresponds to the model (|1|), (|TT]) 
of Phillips, this time specifying the equilibrium variance-density relation in accordance with 
Paveri- Font ana's traffic model. A complete agreement between (^) and ( [Llb| ) results for 
Ce(p, V) = Ce{p), ^e{p, V) = ^e{p), and a special choice of the functional relation r(p, V^)[l — 

p(p,V)]^T{p)[l-p{p)]. 

However, it is not fully justified to assume that the variance Q{r,t) is always in equilib- 
rium Qe{p, V), since the corresponding relaxation time 2/r is of the order of the relaxation 
time 1/r for the velocity V{r,t). Moreover, the approximation 0(r, t) ^ Oe[p(r, t), \^(r, t)] 
does not describe the empirically observed increase of variance B directly before a traffic jam 
develops [0,^. Therefore, we also need the dynamic variance equation (|46|). The remaining 



problem is how to obtain suitable relations for r(r, t) and C(r, t). 

B. Euler-like traffic equations 

Before looking for dynamic relations for T{r,t) and C{r,t), it is plausible first to look 
for equilibrium relations which apply to stationary and spatially homogeneous traffic. For 
this purpose we require the equilibrium solution Pe(f , Vq) of Paveri- Fontana's traffic equation 

Unfortunately, it seems impossible to find an analytical expression for pe{v,vo), but in 
order to derive equations for the velocity moments {v'^) we are mainly interested in, it is 
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sufficient to find the stationary and spatially homogeneous solution pdv) of the reduced 
Paveri-Fontana equation For this we need to know the relation 



Vo{v) = ao + ai5v + 02 {SvY + . . . + a„ )" 



with 



6v := V — V 



(53) 



(54) 



and arbitrary n. However, the equation that determines Vo{v) depends on the unknown 
quantity 



Qo{v) := dvo {vq - Vq) 



2 Pe(^,-t^o) 

Pe{v) 



(55) 



etc. so that we are again confronted with a non-closed hierarchy of equations. 

Luckily, from empirical data and microsimulations we know that the equilibrium velocity- 
distribution 



Pe{v) 



Pe[V) 



(56) 



(at least in the range of stable traffic without stop-and-go waves) is approximately a Gaussian 
distribution |0|,|l],|,||,|3l : 



1 Q-(v-Vef/i2ee) 



(57) 



Inserting (^3]) and (|57D into the equation 



d f ^ Vo{v) — 



dv 



T 



(58) 



which corresponds to equation (|29| ) in the stationary and spatially homogeneous case, we 
find the condition 

d f^Vo{v)-v\ Vo{v)-vdpe , pe f dVo{v) 



dv 

pe 



Pe 



dv ^ T \ dv 



- 1 



[a, - 1) + ( 2a2 - %-}^) 6v + fsaa - {6vf 



■■■ + [kak - ) {6v) 

-(1 -p)pePe 5V. 



Qe 
k~l 



(59a) 
(59b) 
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A comparison of the coefficients of (Sv)^ in ( p9a| ) and (|59b| ) leads to 

a„ = , a„_i = , ... 02 = , cti = 1 , (60) 

and 

ao = K + r(l - p)peee = Vo , (61) 

where we have utilized relation with (|43|) . Consequently, for equilibrium situations 
velocity distribution (|^) implies 

Vo{v) = Vo + 6v. (62) 
With (|57|) and (^2]) we can now derive equilibrium relations for C and F. One obtains 

Te = (63) 

and 

Ce = e,. (64) 

Next, we are looking for relations for non-equilibrium cases. Assuming that the velocity 
distribution 

Pf +\ Pir,v,t) 

P(.;r,t):=-^ (65) 

locally approaches the equilibrium distribution Pe[V{r, t), 9(r, t)] very rapidly, we can apply 
the zeroth-order approximation of local equilibrium: 

P{v- r, t) ^ P(o) {v- r, t) ■= P,[V{r, t), e(r, t)] 

- 1 ^-\v-V{r,t)]y\2e(r,t)] /gg^j 

^27r0(r,t) ■ ^ ^ 

Furthermore, in order to fulfil the compatibility condition 

C(r, t)= [ dv[v - Vir, t)] [vq - Vo{v; r, t)]P{v; r, t) (67) 
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(cf. (|47|)), we must generalize relation ( |62D to 



Voiv;r,t) = Vo + ^^^Sv (68) 



which is fully consistent with (|64D. Relations ( pGf ) and ( pSf ) yield zeroth-order relations for 
the spatio-temporal variation of C(r, t) and J7'(r, t): For the flux density of velocity variance 
we find 

J{r, t) ^ J(o) (p, V, 9) = pr(o) (p, e) = , (69) 

whereas for the covariance the djTiamic equation 

— + = -C- — + - 00 -C -2{l-p)pC\ - (70) 

at or or p or T \ n 



can be derived from the reduced Paveri- Font ana equation (|29D due to 

dv / dfo (5f )^5fop(r, t>, t>o, t) 



dvi6vy[Voiv)-Vo]pir,v,t) 

dv{6vf^p{r,v,t)=J^ (71) 
((5fo := fo — Vo). (The somewhat lengthy but straightforward calculation is presented in 



Ref. [^|. 



In the zeroth-order covariance equation (|70D the quantity 



eo(r, t) := jdvj dvo [vo - Vo{r, t)]'^^^^^^^ (72) 

denotes the variance of desired velocities. The term —QdVo/dr normally vanishes since the 
average desired velocity Vq is approximately constant almost everywhere (cf. ||7^). Due to 
(|6^), the equilibrium variance related to stationary and homogeneous traffic is obviously 
determined by the implicit relation 

ee(Pe, e,) = Ce(pe, Ve, Qe) = Qo ' 2t{1 - p)peQe\ — • (73) 

V TT 

Inserting the above results into equations (^Oj), (0), and (^), we obtain the following 
zeroth-order approximations of the density-, velocity-, and variance-equation respectively: 
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dt dr dr r 

). (76) 

p or T 



Equations (|7iD, (|75D, and ([76| ) are the 'Euler-like equations' of vehicular traffic 



In comparison with the Euler equations for ordinary fluids [p4-|55|| they contain additional 
terms: 

1. The terms (Vq — V)/t and 2(C — 0)/r arise from the acceleration of vehicles towards the 
drivers' desired velocities fo, i.e. they are a consequence of the fact that driver-vehicles 
units are active systems. 

2. The term —(1 —p)pQ results from the vehicles' interactions. It would vanish if mo- 
mentum would be a collisional invariant during vehicular interactions like this is the 
case for atomic collisions |7^. However, without this term the 'vehicular fluid' would 



speed up at bottlenecks which is, of course, unrealistic. 

Moreover, the covariance equation ([70|) is a complementary equation which arises from the 
drivers' tendency to move with their desired velocities vo- 



C. Equilibrium relations and fundamental diagram 

For vehicular traffic, the only dynamic quantity that remains unchanged in a closed sys- 
tem (i.e. a circular road) is the average spatial density p (due to the conservation of the 
number of vehicles). As a consequence, the equilibrium traffic situation is uniquely deter- 
mined by p which obviously agrees with the equilibrium density pe- Equilibrium relations 
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for the average velocity K(pe) and the velocity variance Oe(pe) in dependence of pe = p 
can be obtained from equations (|^) and ( |73D if the relations p{p, V, 9) and r(p, V, 9) are 
given (cf. 0,0]). A simple procedure for finding a solution of these implicit equations is to 

numerically integrate the equations 
dV 

— = Ve[p,,V{y),Q{y)]-V{y) 

= Vo- T{p,, V, 9) [1 - Pip,, V, e)]g,e - V , (77) 
^ = 9e[pe,V(z/),9(y)]-9(i/) 

= 9o - 2r(pe, V, 9)[1 - p(pe, V, e)]ge<S>\ - - 9 (78) 

V 71 

until dV/dy = and dQ/dy = 0. Here, we have replaced pe by = ge{pe,V) in accor- 
dance with section VII B in order to take into account the space requirements of vehicles. 
The theoretical results for the equilibrium velocity- density relation Ve(pe) = lim V{y), the 
equilibrium variance- density relation 9e(pe) = lim Q{y), and the fundamental diagram 

QeiPe) ■= PeVeiPe) (79) 

can be directly compared with empirical data. 

If, however, p(p, V, 9) or r(p, V, 9) are unknown relations, it is still possible to derive 
from the fundamental diagram qdpe) the equilibrium variance-density relation 9e(pe) for 
which an empirical relation seems to be missing: From (|7^) and (|79D we get 

r(l - p)geQe{Pe) = K, " K(Pe) = V, - ^ . (80) 

pe 



Inserting this into ([75| ) we find 

'9e(pe 



9e(Pe) = 9o - 2r(l - p)ge<de{Pe 



9o-2[\/o-K(Pe)]J^^. (81) 



This results in a quadratic equation for the standard deviation yQeiPe) of vehicle velocities 
which is solved by 



a(p.).-^i^,Jfii^M,e„. (82) 

TT V vr 
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V. APPROXIMATE SOLUTION OF PAVERI-FONTANA'S TRAFFIC EQUATION 



The traffic equation of Paveri-Fontana was mathematically investigated in several papers 
dealing with the existence, uniqueness, and numerical determination of a solution which 
satisfies the non-linear initial- value boundary problem [pOf-P^. However, the approximate 
dynamic solution of the reduced Paveri-Fontana equation which will be presented in 
this section has not been proposed before. 

As one would expect, in non-equilibrium situations the zeroth-order approximation (|66D 
does not solve the reduced Paveri-Fontana equation (pOf) exactly. Therefore, we write 



p(r, V, t) =: p(o) (r, v, t) + p(i) (r, v, t) (83) 

with 

P(0)(r,t;,t) := p(r, t)P(o)(.;; r, t) = e-t^-^^'^'*)!^/^®^'-'*)] (84) 

and try to derive a relation for the deviation p(i)(r, f,t). Utilizing that the correction term 
P(i)(r, t>,t) will usually be small compared to p(o)(r, t>,t) we have 

P(i)(r,w,t) < p(o)(r,?;,t) (85) 



and get 



dt ^ dr dv \ T J 
dt dr dv \ T ) 

dt dr T dv T \ dv 



(86) 



(For a detailled discussion of this approximation cf. [^p3t^ . ) Now, introducing the ab- 
breviation 

d d d , ^ 

we can write 
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dp{o) dpjo) ^ ^P(o) ^ ^P(o) dp dp(Q) dV dp(^o) dO 
" ~ ~ dp dt dV dt de dt 



dt dr dt 

_ P(o) dp ^ m ^pm( 



0? 

/ dt 



(88) 



p cit e dt 29 V 
Relations for dp/dt, dV/dt, and dO/dt can be obtained from the Euler-like equations ( [TiD 
(0), and (|7|) via 



We find 



dt dt dr dr 



dp ^ dp dV 
dt dr ^ dr 



(89) 



(90a) 



and 



— -SV—-2Q— -(C-Q) 
dt dr dr r 



(90b) 



(90c) 



For the interaction term we apply a linear approximation in P(i)(r, v, t) which is justified by 
relation (E&). The result is 



(1 — p)p{r, V, t) J dw {w — v)p{r, w, t) 

{1 - p)p(o){r,v,t)p{V - v) - / dw L{v,w;r,t)p{i){r,w,t) 



(91a) 



where we have introduced a linear operator L with the components 



L{v, w; r, t) := (1 - p)p{r, t){[v - V{r, t)]5{v - w) + P(o){v; r, t){v-w)}. 



(91b) 



Here, S{v — w) denotes Dirac's delta function. The linear operator L possesses an infinite 
number of eigenvalues 1/r^ (cf. |]55| , |83| -p^ ) . The relevant eigenvalue is the smallest one since 
it characterizes temporal changes that take place on the time scale we are interested in. It 



is of the order of the average interaction rate per vehicle ||53| , |52| , |54 
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1 _ 1-p 



dv / dw \w — v\p{r,w,t)p{r,v,t) 



w<v 



[l-p)p{r,t)\ — . 

V 71 



The other eigenvalues are somewhat larger 



(91c) 



6| (i.e. < tq for yU 7^ 0) and they 



describe fast fluctuations which can be adiahatically eliminated |78|. As a consequence, we 



can make the so-called 'relaxation time approximation' |87 



dw L(v, w; r, t)p(i)(r, w, t) 



Pii)ir,v,t) 



To 



(91d) 



Now, we calculate 

}Mzl^ + ^ f ^ - 1") - (1 - P)mpiv - V) 



T dv 
1 ( C 

^ 1^0 + 

P(0) 

re 



r 



dv 



+ 



P{o) (C 



T ve 



1) + {^-p)pio)P^v 



(C-0) {l-^^]-{Ve-V)5v 



e 



(92) 



Inserting 
finally obtain 



and (pm to 



into the reduced Paveri-Fontana equation (EDI) we 



P(o) 
P 

P(o) 



dv— 



P- 



dV_' 
dr 



Sv (dv— - _ ^ + l(y _ F) 
G \ dr p dr dr T 



P(o) / (5^^) 
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I \ dr dr T 



-0) 



e 



-P(o)T-o 



\6v 



f 3 6v\ dQ 



1 + 



V,-V 



Sv 



T 



292 



26 / dr 



(93) 



Obviously, the correction term p(i)(r, f,t) is a consequence of the finite interaction free 
time To which causes a delayed adjustment of p{r,v,t) to the local equilibrium p(o)(r, t;, t). 
However, in order to take into account the effects of finite reaction time and braking time we 
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must add a time period r' > to the interaction free time tq. Hence, tq must be replaced 
by the adaptation time 

n^To + r'. (94) 



VI. NAVIER-STOKES-LIKE TRAFFIC EQUATIONS 

With the corrected phase-space density 

p{r, V, t) a; p(o) (r, v, t) + (r, v, t) 

{dvf 3 6v\dQ' 
292 ~ 'dr 

we can calculate corrected relations for the collective ('macroscopic') quantities 



P{Q){r,v,t) 



1 - T* 



where 



We find 



and 



(95) 



F(r,t) ^ {f{v)) := / dvf{v)P^^ ^ F^o){r,t) + F^,)(r,t) (96) 



F(,)(r,t) ^ (/(^))(,) / dvf{v)^^j^ . (97) 



P('^, ^) ~ P(o) (r, t) , V{r, t) ^ V(^o) (r, t) , 
e(r,t) e(o)(r,t) , V{r,t) ^ P(o)(r,t) , (98) 



C{r,t) ^ qo){r,t) . (99) 
However, for the flux density of velocity variance we get 

J{r, t) ^ J^,){p, V, 9) = pr(i)(p, V, 9) = -K^ , (100) 

where 
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K := 3pr,e (101) 



is called a 'kinetic coefficient'. Therefore, the macroscopic traffic equations (|40|) , (p|), and 
(1^) assume form 

dp d(pV) , , 

i + ^ = 0, (102) 

dV dV IdV I 



and 



dt dr p dr pdr \ dr J 

+ -(C-0) + (l-p)«:^. (104) 



Additionally, the corrected covariance equation becomes 

dC ^BC BV VdVo Id f ,de\ 

h V — = -C H C — 

dt dr dr p dr p dr \ dr J 



with the kinetic coefficient 



+ ^ [Ce(p, V, e, c)-c] + ^^^C^ (105) 



c 

(■= K— = 3pr,C (106) 



e 



and the equilibrium covariance 



Ce{p, V, e, C) := Bo - 2r(l - p)pCW - . (107) 



TT 



(For a detailed derivation of ([T0|) to ([TO^) cf. Ref. |7| 



Equations (|102| ), ( |103| ), and ( p,04| ) are the Navier-Stokes-like traffic equations [Q. Com- 
pared with the Navier-Stokes equations for ordinary fluids they possess the additional terms 
{Ve - V)/t and 2{Qe - 0)/r with Be = C + (r/2)(l - p)ndQ/dr which are due to accel- 
eration and interaction processes. Because of the spatial one-dimensionality of the consid- 
ered traffic equations, the velocity equation (|103|) does not include a shear viscosity term 
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{1 / p)d / dr{vQdV / dr) . The variance equation (|104| ) is related to the equation of heat conduc- 
tion. However, G does not have the interpretation of 'heat' but only of velocity variance, 
here. Finally, the Navier-Stokes-like traffic equations are complemented by the additional 
covariance equation ( |105| ) arising from the tendency of drivers to get ahead with a certain 
desired velocity Vq. 

We recognize that the first-order macroscopic traffic equations ( |102| ), (|103|) , ( p.04|) , and 
( |105|) build a closed system of equations. Moreover, according to (p8l), the relations for the 
spatial density, average velocity, velocity variance, and traffic pressure did not change. In 
this sense, the chosen Chapman-Enskog method for closing the hierarchy of macroscopic 
equations is consistent with its assumption, according to which only the expressions for the 
flux density of velocity variance J = pT and the covariance C were to be improved by the 
non-equilibrium correction p(i)(r, However, note that another relation for Vq{v) than 
(|68|) would have led to modiflcations of p, V, and/or 0. 

We also recognize that the flnite adaptation time for approaching the equilibrium 
distribution causes a flnite skewness 

,.= ^ = ^ = = (108) 

' ■ Q3/2 p@3/2 p@3/2 dr y/Q dr ^ ^ 

of the non-equilibrium velocity distribution 

P(.;r,t)^^(°^^^'"'^j+g^)^^'"'^^ (109) 
This leads to the so-called transport terms 

-K— and -CjT- 

or or 

The effect of these terms in equations ( |104| ) and ( |105| ) is to smooth out sudden changes of 
variance and covariance via second spatial derivatives of 0(r, t), namely 

^ (k^^] and ^ (c^^] (111) 
dr \ dr I dr \ dr I ' 



30 



VII. CORRECTIONS OF THE MODEL 



A. Driver behavior and bulk viscosity 

We remember that the term —{l/p)dV/dr describes an anticipation effect. It reffects 
that drivers accelerate when the 'traffic pressure' V = pO lessens, i.e. when the density p or 
the variance G decreases. However, drivers additionally react to a spatial change of average 
velocity. This effect can be modelled by the modified pressure relation 

dV 

V{p,V,Q):=pQ-v-Q^ (112) 



which gives velocity equation ( |103| ) a similar form like variance equation ( |104D and covariance 
equation (|105|) . 



In order to present reasons for relation (|112|) let us assume that drivers switch between 



two driving modes m G {1,2} depending on the traffic situation. Let m = 1 characterize 
a brisk, m = 2 describe a careful driving mode. Then, we can split the density p(r, t) into 
partial densities pm{r,t) that delineate drivers who are in state m: 

pi{r,t) + p2{r,t) = p{r,t) . (113) 

Both densities are governed by a continuity equation, but this time we have transitions 
between the two driving modes with a rate R{pi, V) so that 

^ = -|:(Pi^)-^(Pi.^), (114a) 
^ = -^iP2V)+Rip-p2,V). (114b) 

Adding both equations we see that the original continuity equation (|102|) is still valid. Now, 
defining the substantial time derivative 

--# + ^# (118) 

Dt dt Or ^ ^ 

we can rewrite ( |114a| ) and obtain 
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D / Dt describes temporal changes in a coordinate system that moves with velocity V . As- 
suming that pi relaxes rapidly we can apply the adiahatic approximation ||7^ 

Dpi 



Dt 







(117) 



which is valid on the slow time-scale of the macroscopic changes of traffic flow. This leads 
to 

dV 



dr 



(118) 



Relation (|1 1 7|) implies that the density pi of briskly behaving drivers is approximately con- 
stant in the moving coordinate system whereas the density p2 = p — pi of carefully behaving 
drivers varies with the traffic situation: 



DP2 

Dt 



dV 



(119) 



p2 increases when the average velocity spatially decreases {dVjdr < 0) since this may 
indicate a critical traffic situation. 

According to relations (|1 14|) , (|118|) incessant transitions between the two driving modes 
take place as long as traffic flow is spatially non- homogeneous (i.e. dV/dr ^ 0). This leads 
to corrections of the pressure relation. Expanding V with respect to the variable R which 



characterizes the disequilibrium between the two driving modes we find |74 



p(p,e,i?) = p(p,e,o)- — 



R=0 



dV_ 

P^ Qj. 



(120) 



With the equilibrium relation V{p, 6, 0) = p6 and 

dV 



V ■= Pi 



dR 



;i2i) 



R=0 



we finally obtain the desired result 



Vip,Q,R)=Vip,V,Q) = pQ-ri 



dV 

dr 



(122) 



A more detailed discussion can be found in Ref. [74 
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B. Modifications due to finite space requirements 



We will now introduce some corrections that are due to the fact that vehicles are no 
point-like objects but need, on average, a space of 

s{V)=l + VT (123) 

each. Here, I > Iq is about the average vehicle length whereas VT corresponds to the safe 
distance each driver should keep to the next vehicle ahead. T is about the reaction time. 
Consequently, if AN{r, t) :— p{r, t) Ar means the number of vehicles that are at a place 
between r — Ar/2 and r + Ar/2, the effective density is 

t) = ^^(^'^^ = (124) 

' ' l^r-m{r,t)s[V{r,t)] I - p{r,t)s[V{r,t)Y ^ ' 

Since AA^(r, t)s{V) is the space which is occupied by AN{r, t) vehicles, the effective density 
is the number AA'"(r, t) of vehicles per effective free space Ar — AA^(r, t)s{y). 

The reduction of available space by the vehicles leads to an increase of their interaction 
rate. Therefore, we have 

■= (^~P) J J \v - w\Q{r, V, Wo, t)p{r, w, vq, t) 

V / tr y 

V 

— (1— p) j dw j dwo\w — v\g{r,w,wo,t)p{r,v,vo,t) (125) 



with 

g{r, V, vo, t) := . 126 

1- p{r,t)s[V{r,t)] 

Consequently, we obtain the corrected relation 

i:=(l-p)m/®. (127) 

To V TT 

In addition, we must replace V and J' by 
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respectively ||76[. For the kinetic coefficients t], k, and ( we obtain the corrected relations 



and C' := z ^ = 3f?r,C . (129) 

l-ps{V) 



The corrected formula 



for the equilibrium pressure corresponds to the pressure relation of van der Waals for a 'real 
gas'. According to ( |130| ), the traffic pressure diverges for p p^ax '-=1/1 which causes a 



deceleration of vehicles. 

The corrected kinetic coefficients ri'{p,V,Q), k'{p,V,Q), and ('{p,V,Q,C) also diverge 
for p pinax We find for example 



n 3^re = ^-^ (131) 

so that the divergence of n' is a consequence of the finite reaction- and braking-time r'. 
This divergence causes a homogenization of traffic flow since the second spatial derivatives 
d/dr{ridV/dr), d/dr{K,dQ/dr), and d/dr{(dQ/dr) produce a spatial smoothing of average 
velocity V, variance B, and covariance C respectively. 

It is the divergence of 'traffic pressure' and kinetic coefficients for p Pmax that prevents 
the spatial density p from exceeding the maximum density pmax @1- 



VIII. SUMMARY AND OUTLOOK 

This paper started with a discussion of the most widespread macroscopic traffic models. 
Each of them is suitable for the description of certain traffic situations on freeways but fails 
for others. Therefore, an improved fluid-dynamic model was derived from the gas-kinetic 
traffic equation of Paveri- Font ana which is very well justifled and does not show the 
peculiar properties of Prigogine's Boltzmann-like approach [Q . 
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For the derivation of the improved traffic model, moment equations for collective ('macro- 
scopic') quantities like the spatial density, average velocity, and velocity variance had to be 
calculated. The system of macroscopic equations turned out to be non-closed so that a 
suitable approximation was necessary. Here, the well proved Chapman-Enskog method was 
applied. In zeroth-order approximation the velocity distribution is assumed to be in 'local 
equilibrium'. According to empirical data, the latter is characterized by a Gaussian velocity 
distribution. Depending on the respective kind of zeroth-order approximation one arrives 
at the Lighthill-Whitham model |]10|, the model of Phillips [p|,[7l|, or the Euler-like traffic 



equations. 

For the derivation of a first-order approximation, the reduced Paveri-Fontana equation 
was linearized around the local equilibrium solution and solved by application of the Euler- 
like traffic equations. The resulting correction term for the non-equilibrium velocity distri- 
bution allowed the calculation of additional transport terms which describe a flux density of 
velocity variance and covariance in spatially non-homogeneous situations. They are related 
with a finite skewness of the velocity distribution. The shear-viscosity term vanishes because 
of the one- dimensionality of the considered traffic equations. Nevertheless, a bulk-viscosity 
term results from transitions between two different driving modes: a brisk and a careful one. 

The resulting Navier-Stokes-like traffic equations were finally corrected in order to take 
into account the finite space requirements of vehicles. They overcome the deficiencies of the 



former macroscopic traffic models so that the criticism by Daganzo and others could be 
invalidated: 

1. The anticipation term which, in other models, is responsible for the prediction of 
negative velocities vanishes in problematic situations like the one described at the end 
of Sec. II since the variance becomes zero, then. 

2. The density p{r,t) does not exceed the maximum admissible density pbb (= bumper- 
to-bumper density) [^] since the diverging viscosity term causes a homogenization 
of traffic flow and the diverging traffic pressure suppresses an unrealistic growth of 
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velocity which stops a further increase of traffic density. 

3. The model takes into account different driving styles by a distribution of desired ve- 
locities Vo which are directly associated with the individual drivers. An extension of 
the Navier-Stokes-like traffic equations to different vehicle types (cars and trucks) is 
possible [g. 

4. The interaction between drivers is modelled anisotropically since the slower vehicle is 
assumed not to be affected by a faster vehicle behind it or overtaking it. 

5. According to the Navier-Stokes-like equations, disturbances may propagate with a 
velocity c > V since a certain proportion of vehicles moves faster than the average 
velocity V due to the finite velocity variance B. Therefore, in contrast to what was 
claimed by Daganzo it is admissible that macroscopic traffic models "exhibit one 



characteristic speed greater than the macroscopic fluid velocity" ||5T|j89|] . 



Present investigations focus on the computer simulation of the Navier-Stokes-like traffic 
equations. This work has already been successfully started for a circular road PPPj] and is 
now extended to complex freeway networks. 

Moreover, the gas-kinetic and Navier-Stokes-like traffic models can be generalized to 
models for multi-lane traffic where overtaking and lane-changing is explicitly taken into 



account [^. By this, formulas for the relations r(p, V, 9) and p{p, V, 9) can be derived pT 
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